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ABSTRACT 

We report on INTEGRAL observations of the bright black-hole transient GX 339—4 
performed during the period August-September 2004, including the fast transi- 
tion (10 hr) observed simultaneously with INTEGRAL and RXTE on August 15th 
and previously reported. Our data cover three different spectral states, namely 
Hard/Intermediate State, Soft/Intermediate State and High/Soft State. We investi- 
gate the spectral variability of the source across the different spectral states. The hard 
X-ray spectrum becomes softer during the HIMS-to-SIMS transition, but it hardens 
when reaching the HSS state. A principal component analysis demonstrates that most 
of the variability occurs through two independent modes: a pivoting of the spectrum 
around 6 keV (responsible for 75% of the variance) and an intensity variation of the 
hard component (responsible for 21%). The pivoting is interpreted as due to changes 
in the soft cooling photon flux entering the corona, the second mode as fluctuations 
of the heating rate in the corona. These results are very similar to those previously 
obtained for Cygnus X-1. Our spectral analysis of the spectra of GX 339—4 shows a 
high energy excess with respect to pure thermal Comptonisation models in the HIMS: 
a non-thermal power-law component seems to be requested by data. In all spectral 
states joint IBIS, SPI and JEM-X data are well represented by hybrid thermal/non- 
thermal Comptonisation (eqpair) . These fits allow us to track the evolution of each 
spectral components during the spectral transition. The spectral evolution seems to 
be predominantly driven by a reduction of the ratio of the electron heating rate to the 
soft cooling photon flux in the corona, Ih/ls- The inferred accretion disc soft thermal 
emission increases by about two orders of magnitude, while the Comptonised luminos- 
ity decreases by at most a factor of 3. This confirms that the softening we observed 
is due to a major increase in the flux of soft cooling photons in the corona associated 
with a modest reduction of the electron heating rate. 
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1 INTRODUCTION 

Black hole candidates (BHCs) are known to show different 
spectral states in the X and 7-ray domain. The two main 
spectral states are the Low/Hard State (LHS), with the high 
energy spectrum described by a cut-off power-law (typically 
r ~1.5 and Ecut ~100 keV), and the High Soft State (HSS) 



* Based on observations with INTEGRAL, an ESA project with 
instruments and science data centre funded by ESA member 
states (especially the PI countries: Denmark, France, Germany, 
Raly, Switzerland, Spain), Czech Republic and Poland, and with 
participation of Russia and the USA. 



with a thermal component peaking at few keV and the hig h 
energy power-law much softer (F > 2.2) I Zdziarski 2000l ). 
Usually, this spectral variability is interpreted as due to 
changes in the geometry of the central parts of the accre- 
tion flow (see Done, Gierliriski & Kubota 2007 for a recent 
review) . In the L HS the standard geometrica lly thin and op- 
tically thick disc ( Shakura fc Sunvaev 19731 ) would be trun- 
cated far away from the last stable orbit. In the innermost 
parts, a hot accretion flow is responsible for the high energy 
emission, via thermal Comptonisation of the soft photons 
coming from the truncated disc. In the HSS, the optically 
thick disc would extend close to the minimum stable orbit 
producing the dominant thermal component. The weak non- 
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thermal emission at higher energy is believed to be due to 
up-scattering of the soft the rmal disc emission in active co ro- 
nal regions above the disc (jZdziarski fc Gierkihski 2 004). 

However, recently long XMM-Newton observations of 
GX 339—4 have shown that a standard thin accretion disc 
may remain at the innermost stable circular orbit around 
the black hole also during the LHS (Miller et al. 2006, Reis 
et al. 2008) 

Since its discovery IjMarkert et al. 19731 ). the X-ray bi- 
nary GX 339—4 has been thoroughly studied at all wave- 
length from radio to gamma-rays. Its star companion is still 
unknown, even though upper limit on the optical luminos- 
ity allowed to classify the source as Low Mass X-ray Binary 
(LMX B; Shahbaz, Fender & Charles 2001). Classified as 
BHC (jZdziarski et al. 19981 ) . GX 339—4 is a transient source 
spending long periods in outburst. Hynes et al. (2003) esti- 
mated a mass function of 5.8 ± 0.5 Mq and Zdziarski et al. 
(2004) derived a lower limit on the distance of 7 kpc. Before 
the launch of RXTE the source had been observed mostly 
in the Low/Hard spectral state (LHS), even though few 
spectral trans itions to the Very Hi gh state (VHS) were re- 
ported before (JMivamoto et al. 1991 ). Thereafter GX 339-4 
remained bright and mostly in the LHS until 1999 when it 
went into quiescence. After the quiescent state observed by 
BeppoSAX (Kong et al. 2000, Corbel et al. 2003), GX 339-4 
showed two new outbursts: in 2002/2003 (Smith et al. 2002a, 
Belloni et al. 2005) and in 2004 after one year in quiescence 
(Buxton et al. 2004, Belloni et al. 2004). The long-term vari- 
ability of GX 339-4 (1987-2004) is extensively presented in 
Zdziarski et al. (2004). 

The spectral evolution of the outburst of a black hole in 
low-mass system is often investigated by plotting the flux of 
the source as a function of the X-ray hardness. In this Hard- 
ness Intensity Diagram (HID), GX 339—4 usually follows a 
q-\ike pattern throughout the outburst (Homan & Belloni 
2005; Belloni 2005). The HSS and LHS correspond respec- 
tively to the left and right hand side vertical branches of the 
q. State transitions occur when the source crosses the upper 
or lower horizontal branches, those constitute intermediate 
states between the LHS and HSS. Depending on the location 
of the source on a horizontal branch, Belloni et al. (2005) 
have identified two different flavours of intermediate state: 
the Hard Intermediate State (HIMS) and Soft Intermediate 
State (SIMS). An alternative states classification defines the 
upper horizontal branch as VHS (see Tanaka & Lewin 1995; 
van der Klis 1995; McChntock & Remillard 2006). 

During a typical outburst the source starts in a faint 
LHS, move upward along the LHS branch, then move left- 
ward towards the HSS along the VHS branch, and when 
luminosity decreases the source moves down along the HSS 
before transiting right back to the LHS along the lower hor- 
izontal branch, and then back to quiescence in the LHS. 
A puzzling property of this evolution is that the spec- 
tral changes in GX 339—4 lag the variations of the lu- 
minosity: hard-to-soft state transitions during the rising 
phase occur at higher luminosities than the soft-to-hard 
ones during the d eclining phase. The so-called "hystere- 
sis" IjZdziarski fc G icrkinski 2004) is also observed in other 
LMXBs such as GRS 1758-258 and IE 1740.7-2942 (Smith 
et al. 2002b; Del Santo et al. 2005). 

The aim of this work is to use INTEGRAL 
(jWinkler et al. 20031 ) data collected during the softening of 



the 2004 outburst to study the broad-band spectral evo- 
lution of GX 339—4 during a hard-to-soft state transition. 
Some of these data obtained on 2004 August 14*''-16"', right 
across the HIMS-SIMS fast transition, are presented in Bel- 
loni et al. (2006) . This transition, marked by the disappear- 
ing of type-C QPO and by the appearance of a type-B QPO 
(Nespoli et al. 2003), is also associated to the ejection of fast 
relativistic jets, which has led to the identification of a 'jet- 
line' in the HID (Fender, Belloni & Gallo 2004). Moreover, 
focussing on the high energy part of the spectrum, Belloni 
et al. (2006) reported on the disappearing of the high energy 
cut-off (measured at 70 keV in the HIMS) in the SIMS. Re- 
cently, a similar evolution of such cut-off has been found for 
the BHC GRO J1655-40 (jjoinet et al. 2008. ). These authors 
propose that the cut-off increased significantly or vanished 
completely as the radio-jet turned off in the SIMS. 

Here, we extend the work preliminarily presented in Del 
Santo et al. (2006) by studying the spectral evolution of 
GX 339—4 over a longer period of time covering almost the 
whole span of the LHS-HSS transition (from 2004 August 
9 to 2004 September 11). In this work, different Compton- 
isation models for spectral fitting have been used; further- 
more, a principal component analysis provided the variabil- 
ity modes of GX 339-4. 



2 OBSERVATIONS AND DATA ANALYSIS 

The 2004 X-ray activity of GX 339-4 started on Febru- 
ary 9"" (Smith et al. 2004, Belloni et al. 2004). This out- 
burst was observed two times by the INTEGRAL satel- 
lite: during the first part of the hard X-ray activity (19 
February-20 March) and during the decay of the hard X- 
ray fiux, i. e. 9 August-11 September (the evolution of the 
fiux and hardness of the source during this second period 
is shown in Fig[T]). As expected in black hole transients 
(jHoman fc Belloni 20051 ). at the beginning of the outburst 
GX 339-4 was in the Low/Hard state {INTEGRAL obser- 
vations presented in Joinet et al. 2007). On 2004 August 
15 {INTEGRAL orbit 224), a fast spectral state transition 
Hard-Intermediate (HIMS) to Soft-Intermediate (SIMS) was 
reported u sing simultaneous RXTE and INTEGRAL ob- 
servations (|Belloni et al. 20061 ). Thereafter GX 339-4 dis- 
played a softer state as also shown by the hardness ratio 
(Figdl bottom panel). 

This paper is based on INTEGRAL public data col- 
lected during the second part of the outburst (including 
the fast spectral state transition) between revolution 222 
and 233. Because of the smaller JEM-X FOV (4° x 4°) 
with respect to IBIS (9° x 9°, fully coded), we selected 
a data sub-set cont aining a ll common pointings of JEM- 
X ("Lund ct al. 20031) IBIS (|Ubertini et al. 2003h and SPI 
(Vcdrcnnc ct a l. 20031 ). for a total of 76 Science Windows 
(SCWfEl) (see Table [7] for details). 

We reduced data of the IBIS low energy detector, IS- 
GRI l|Lebrun et al. 20031 ). and JEM-X data using the IN- 
TEGRAL Off-Line Scientific Analysis (OSA; Courvoisier et 
al. 2003). 

The SPI data analysis has been performed using a 



^ Each sew corresponds to a pointing of about 30 minutes. 
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Figure 1. From top to bottom: GX 339—4 counts rate in 1.5-3 
keV, 3-6 keV, 20-40 keV with ASM, JEM-X and IBIS/ISGRI, re- 
spectively. The related INTEGRAL orbits are also marked on the 
IBIS light curve. IBIS/ISGRI to JEM-X hardness ratio is shown 
in the bottom panel. The dotted line marks the fast transition 
HIMS-to-SIMS caught by Belloni et al. (2006). 



model fitting algorithm (JBouchet et al. 2 003). With the 
sources emitting in hard X-ray present within about 15 de- 
grees around GX 339—4, we built a sky model consisting 
of 6 sources (4U 1700-377, OAO 1657-415, 4U 1630-47, H 
1636-536, IGR J16318-4848 and IGR J16320-4751), consid- 
ered as constant on over a revolution, except 4U 1700-377, 
for which a ~2000 seconds variability time-scale was used. 
The background pattern was built from the revolution 220. 

The broad-band spectra were averaged in seven groups 
corresponding to the periods in Tab. [T] A division of the 
different groups into spectral states based on timing analysis 
of RXTE data, i. e. on the type of the observed QPOs, is 
also presented in the same table. 

Because of the low statistics, SPI averaged spectra have 
been used only for spectral fitting of periods 1, 2 and 3. 
JEM-Xl camera spectra have been averaged in the energy 
range 4-20 keV for periods 1, 2, 3 and 4, and in the band 
3-13 keV for 5 and 7. ISGRI and SPI data above 20 keV and 
22 keV, respectively, have been used, while the high energy 
limits were chosen between 200 and 600 keV depending on 
the statistics. 

Spectral fits were performed using the spectral X-ray 
analysis package XSPEC v. 11.3.1. 



3 SPECTRAL VARIABILITY 

As a first step in quantifying the spectral evolution of the 
source, we attempted to fit each of the 76 IBIS/ISGRI 
pointing spectra with a simple power-law (powerlaw under 
XSPEC) in the energy band 20-80 keV. Due to the consis- 
tent flux decreasing of GX 339—4 (Fig. [T] panel 3 from the 
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Figure 2. GX 339—4 photon index evolution. The IBIS spectra 
collected sew by sew of periods 1 and 2 have been fitted with a 
simple pawer-law in the range 20-80 keV (triangles). Starting from 
period 3, all spectra have been averaged for each period and fitted 
with same model (diamonds). The HIMS-to-SIMS transition is 
marked (dotted line). 



top), from period 3 onward the statistics by pointing is too 
poor to provide meaningful fit parameters. This is why from 
period 3 to 7, we fitted the spectra averaged over each data- 
set. The evolution of the photon index (Fig. [2]) indicates 
that the spectrum softens when reaching the transition (F 
goes from 1.9 to 2.4). Thereafter it becomes harder when 
the soft state is reached (F ~ 2 in rev 226). The hard X- 
ray spectrum therefore appears softer around the end of the 
transition (i.e. in the SIMS) than in the HSS. 

Then, we fitted the seven averaged broad-band spec- 
tra (IBIS, JEM-X and SPI when available) with a simple 
model combining a multicolor disc black-body (diskbb), 
a pawer-law (with a high energy cut-off when requested) 
and an emission line (gauss) fixed at 6.4 keV. Since the 
estimated absorption for GX 339—4 is rather low and we 
did n ot have low-energy coverage, we fixed it to 5x10^^^ 
cm~^ iMehdez fc van der Klis 1997^ . Best-fit parameters of 
six periods (period 6 was avoided because of the statistic) 
are summarised in Tab. [2] 

The main spectral changes are the steepening of the 
power-law component and disappearance of the high energy 
cut-off during the transition (periods 3 and 4). This is in 
agreement with the results obtained using simultaneous IN- 
TEGRAL and 7?Xri? observations and extensively discussed 
in Belloni et al. (2006). During periods 5 and 7, the spec- 
tral parameters are consistent with being constant and the 
power- law is harder than in the SIMS (period 4). The disc 
temperature was frozen at 0.5 keV for periods 1 and 2, while 
for the others it was possible to constrain the parameter 
kTii. 



4 AVERAGED SPECTRA 

Since the cut-off power-law is simply an empirical model de- 
scribing Comptonised spectra, we have used a further model 
which provides a physical description of X/7-ray spectra. In 
this section, we present spectral fits obtained using a hy- 
brid the rmal/non-th ermal Comptonisation model, namely 
EQPAIR (|Coppi I999I ). and relativistic iron line line emission 
and disc reflection. The best-flt parameters of six periods 
are shown in Tab. [3] 

In the EQPAIR model, emission of the disc/corona sys- 
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Figure 3. Joint JEM-X, IBIS and SPI (only data-sets 1, 2 and 
3 for the latter) energy spectra of GX 339—4 during periods 1 
(black), 2 (blue), 3 (red), 4 (yellow), 5 (green), 7 (orange). The 
data are fitted with the thermal/non-thermal hybrid Comptoni- 
sation model EQPAIR plus DISKLINE (see Tab. [J}- In order to give 
a comparison with a pure Low/Hard state spectrum, we show the 
averaged spectrum collected during rev 175 (violet) and reported 
in Joinet et al. (2007). 



tern is modeled by a spherical hot plasma cloud with con- 
tinuous acceleration of electrons illuminated by soft pho- 
tons emitted by the accretion disc. At high energies the dis- 
tribution of electrons is non-thermal, but at low energies 
a Maxwellian distribution with temperature kT^ is estab- 
lished. 

The properties of the plasma depend on its compact- 
ness / = Lo-t/Rtii^c' where L is the power of the source 
supplied by different components, R the radius of the sphere 
(10^ cm) and or is the Thomson cross-section. Is, kh, Inth 
and Ih — Uh + Inth correspond to the power in soft disc pho- 
tons entering the plasma, thermal electron heating, electron 
acceleration and total power supplied to the plasma. 

The spectral shape is insensitive to the exact value of 
the compactness but it depends strongly on the the compact- 
ness ratios Ih/h and Inth/lh- It is customary, when fitting 
data with this model, to fix Is to some reasonable value and 
parametrise the model with these compactness ratios. Leav- 
ing Is as free parameter we do not obtain any improvement 
on the x^- We therefore fixed Is — 10 following Gierlinski et 
al. (1999). On the contrary, we cannot fix Ih because it is 
not included as free parameter in the eqpair version imple- 
mented in XSPEC. 

If we knew with some accuracy the distance and size 
of the Comptonising region, the normalisation of the spec- 
trum could provide a constraint on the absolute value of 
the compactness. Since, however, there are large uncertain- 
ties on the size of the emitting region, and since there could 
be several active regions in the accretion fiow contributing 
to the observed fiux, it is not straightforward to unambigu- 
ously constrain compactness from the observed luminosity. 
For this reason, in the version of eqpair used to fit the data 
there is an additional normalisation parameter which actu- 
ally disconnects Is (as a fit parameter) from the observed 
luminosity. 



The energy balance in the Comptonising medium de- 
pends mainly on the ratio Ih/h ■ For larger values of Ih/h the 
energy of the electrons is on average higher and, as a con- 
sequence, spectrum from Comptonisation is harder. From 
the best-fit parameters shown in Tab.O we indeed measure 
lower values of Ih/h as the spectra become softer. Since in 
our spectral fits we chose to fix Is = 10, in this model the 
variations of Ih/h are only due to changes in Ih- 

However, due to the free normalisation parameter of 
EQPAIR we do not know if this what really happens in 
GX 339—4 (i.e . h is constant while Ih changes). We know 
for sure that we observe a change in Ih/ls but we cannot 
tell from the fit parameters whether this evolution is due 
to changes in the heating rate of the corona, changes in the 
luminosity of the disc, or both. 

Nevertheless, from the best fitting model we can get im- 
portant insights by computing the absolute disc and Comp- 
tonised fiuxes (respectively Fbb and Fcompt)- If we consider 
the absolute values of Fcompt and Fbi, (see Tab. |3}, we see 
that the thermal disc fiux changes by more than one order of 
magnitude while the Comptonised fiux decreases by a factor 
of about three. The softening we observe is therefore caused 
by a dramatic increase in the disc thermal fiux in the corona 
associated with a modest reduction of the electron heating 
rate. 

The size of the emitting region could be constant, or 
not. Assuming Is constant, as we did, it would imply that 
the size of the emitting regiorjj increases by a factor of ~ 14. 
On the contrary, the timing features usually indicate that 
the emitting region becomes smaller when a source evolves 
from the LHS toward the HSS (see e.g. Gilfanov, Chura- 
zov, Revnivtsev 1999). If this is the case in our observation, 
the soft compactness actually increases during the transi- 
tion. The bolometric luminosity increased by a factor of two 
during the whole spectral evolution. 

The seed photons temperature kT^t was frozen at 300 
eV for the spectra of revolutions 222 and 223; thereafter it 
was possible to constrain this parameter and observe a slight 
increase (Tab. [3)1 simultaneously with the flux of soft cooling. 
This suggests that the increase of soft cooling photons is due 
to the higher temperature of the disc. 

However, the disc luminosity tends to decrease for larger 
temperatures. The highest disc temperatures, achieved dur- 
ing period 5 and 7, do not correspond with the strongest 
observed blackbody fiuxes of period 3 and 4. This suggests 
that the disc emitting area and geometry of the corona are 
not constant during the transition. The radiating disc sur- 
face area appears smaller in the HSS (periods 5 and 7) than 
in the intermediate states (periods 3 and 4). This suggests 
that the corona becomes more compact as the disc tempera- 
ture increases. This would be consistent with the truncation 
radius of the accretion disc moving closer to the black hole 
in the HSS. 

In order to consider the case of a hybrid plasma, the 
model allows for a fraction l„th/lh of the power to be injected 
in the form of non-thermal electrons rather than heating of 
the thermalised distribution. The non-thermal electrons are 
injected with a power-law distribution 'y~''«"j , with Lorentz 



^ or more exactly the number of X-ray emitting regions times 
their typical size 
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factors ranging from 7min = 1-3 to 'jmax = 1000. From our 
spectral fits, we infer dnj values in the range 2-3, as ex- 
pected form shock acceleration models. We found large non- 
thermal fractions in all spectra, including the HIMS ones. 
This is in agreement with previous eqpair fits of Cygnus 
X-1 in intermediate states (Gierlinski et al. 1999; Malzac et 
al. 2006). 

However, we tried to fix Inth/lh = for all spectra. 
It resulted that only for the period 1 we obtained a good 
xt (i.e. 0.98) with remaining parameters consistent within 
the errors to the ones in Tab. [3] For the other spectra, not 
negligible Inth/lh fraction is required. Unfortunately the un- 
certainties on Inth/lh are quite large and do not allow us 
to draw any conclusion regarding its evolution during the 
spectral transition. 

Another interesting free parameter is the Thomson op- 
tical depth (tes). It is related to the ionization electrons 
only. The total optical depth (jt) is the sum of the opti- 
cal depth of e^e~ pairs plus Tea- I n agreement with pre- 
vious studies (|Gierliriski et al. 19991 ) our fits show that the 
electron-positron pairs constitute only a small, if not neg- 
ligible, fraction of the Comptonising leptons. An interest- 
ing trend is that the Thomson optical depth of the corona 
decreases significantly during the transition. This could be 
linked to the 'compactification' of the corona inferred from 
the evolution of the soft photon compactness and disc lu- 
minosity, if the electron density remains constant while the 
size of the corona decreases. 

The disc reflection component is calculated for neutral 
material with standard abundances. We assumed an incli- 
nation angle of 50 deg. This component is then convolved 
with a general relativistic kernel assuming reflection on an 
accretion disc extending from 6 to 1000 gravitational radii 
around a Schwarzschild black hole. 

We checked the need for the reflection component freez- 
ing Q./2-K at 0. The results show that the x a^re consistently 
worst (F-test probabilities < 10~*), confirming that intro- 
ducing the reflection component improves significantly the 
fits. We can conclude that Compton reflection is requested 
by the data. Our fits suggest that the reflection component 
increases up to f2/27r=l as the spectrum goes in softer states. 
This is expected in the truncated disc model as the system 
evolves from a geometry where the reflecting disc is trun- 
cated at a large distance from the black hole (and therefore 
intercepts a small fraction of the coronal radiation) to a 
situation where the accretion disc is sandwiched by the illu- 
minating corona. In period, 4, 5 and 7 the statistics was too 
low to constrain the reflection parameter and we decided to 
freeze it to _R = 1. 

We used the XSPE C model diskline to model the iron 
line (JFabian et al. 19891 ). Due to the energy resolution and 
line sensitivity, the iron line study cannot be performed with 
JEM-X. The line energy was therefore imposed at 6.4 keV. 
The inclination angle, inner radius of the disc and the disc 
emissivity law were fixed at the same values as in the Comp- 
ton reflection model. The normalisation is therefore the only 
free parameter for the iron line emission. 

Spectra and models of the six different periods are 
shown in Fig. [S] For a comparison between all the st ates, 
the LHS spectrum observed in July (jjoinet et al. 2007^ has 
been plotted in the same figure. 
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Figure 4. Residuals obtained with a simple thermal Compton- 
isation model (COMPPS) by fitting IBIS and SPI spectra of rev 
223. 
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Figure 5. IBIS and SPI energy spectra and residuals of rev 223 
fitted with a Comptonisation model (cOMPPS) plus power-law. 



5 THE HIGH ENERGY EXCESS 

During the Low/Hard state of the GX 339—4 outburst oc- 
curred in 1991, a high energy excess (above 200 keV) was 
observed with OSSE (Johnson et al. 1993; Wardziiiski et 
al. 2002). In March 2004, i. e. during the increasing phase 
of the 2004 outburst, Joinet et al (2007) observed a sim- 
ilar feature with SPI when GX 339—4 was in its canoni- 
cal LH state. Similar features have been recently reported 
for other black hole transients in LH states observed dur- 
ing the decreasing phase of the outbursts (Kalemci et al. 
2006; Kalemci et al. 2005). It is also more c ommonly asso- 
ciated with state transitions ([Malzac et al. 2006) and with 
soft states (jGierliiiski et al. 19991 ) ( in which case the non- 
thermal emission actually dominates over thermal Comp- 
tonisation) . 

Such a hard tail is usually attributed to the presence of 
a small fraction of non-thermal electrons in the hot Comp- 
tonising plasma. In the LHS and HIMS such a high energy 
excess could also be the result of spatial/temporal varia- 
tions in the plasma parameters (see e.g. Malzac & Jourdain 
2000). In the HIMS of our observations, i.e. rev 222, 223 
and 224.1, GX 339—4 shows a highly significant evidence of 
a cut-off around 70-100 keV. This indicates that soft photons 
are Comptonised by a thermal electrons population with a 
temperature kT^ that can be measured by the cut-off en- 
ergy. In addition, we found non- negligible values of Inth/lh 
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Energy. keV 



Energy. keV 



Figure 6. PCI {left) and PC2 (right). From top to bottom panels: (i) the shape of the principal component (which would be added or 
subtracted to the average spectrum in order to reproduce the variability); (ii) the effects on the shape and normalisation of the spectrum 
(time-averaged spectrum, shown by the light line, and spectra obtained for the maximum and minimum amplitude of the component); 
(iii) the ratios of the maximum (and minimum) spectrum to average; (iv) the contribution of the component to the total sample variance 
as a function of energy. 



for revolutions 223 and 224.1, indicating that a non-thermal 
emission is also requested by the data. 

In order to investigate furtherly the possible presence 
of a non-thermal component during the HIMS, IBIS and 
SPI spectra of the three HIMS periods (rev 222, 223 and 
224.1) have been fitted with the COMPPS model of Poutanen 
& Svensson (1996) (see Tab.lH). Using this thermal Comp- 
tonisation model, some residuals are present at high energy, 
especially in the spectrum of rev 223 (Fig. |4]). In order to 
mimic the presence of a non-thermal component, we added a 
power-law to the pure thermal Comptonisation model (Fig. 
[5]). This led to an improvement of the y^ that is highly signif- 
icant for orbit 223 (Tab.©. The F-test probability that this 
improvement was by chance is Pcom-pps ~ 6.6 x 10~^. On the 
contrary, because of the low statistics in rev 224.1, adding a 
power-law to the Comptonisation model is not required. 

COMPPS is a highly accurate iterative scattering Comp- 
tonisation model in which subsequent photon scatterings 
are directly followed. In practice this method is limited to 
Comptonisation in a plasma of reasonably small Thomson 
optical depth (r < 3 for spherical geometry). In order to be 
sure that the presence of a non-thermal component does not 
depend on the specific Comptonisation model, we also fitted 
the data with the comptt model of Titarchuk (1994) which 
is based on an approximate solution of the kinetic equation 
with some relativistic corrections. We performed the same 
fits as with COMPPS. As expected, we found significantly dif- 
ferent numerical values of the best-fit parameters (see Tab. 
|4|. However, we obtained a stronger indication regarding 
the presence of the non-thermal component: the addition of 



a power-law leads to a highly significant improvement of the 
X^ both in spectra of rev 222 and rev 223. 



6 THE PRINCIPAL COMPONENT ANALYSIS 

In order to study the overall variability of GX 339—4 we 
present in this section a further analysis, i. e. Principal Com- 
ponent Analysis (PC A), based on a different approach. The 
PCA is a powerful tool for multivariate data analysis used 
for a broad range of appli cations in natural as well as so- 
cial science IjKendall 198Ci ). The main use of the PCA is 
to reduce the dimensionality of a data-set while keeping as 
much informations as possible. PCA transforms a number 
of (possibly) correlated variable in a smaller number of un- 
correlated variables called principal components. 

Details on the method, applied to the analysis of the 
spectral variability of Cyg X-1 with INTEGRAL, are given 
in Malzac et al. (2006). It provides us an approximation of 
the broad-band energy spectrum for each pointing that will 
be used for a physical interpretation of the variability. This 
method is much more convenient than fitting the spectra 
for each SCW, also considering the large number of model 
parameters and the poor statistic of spectra in short (3-4 
ks) exposure times. 

Here, we applied a similar procedure to the state transi- 
tion data of GX 339—4. We used p=7& spectra (one for each 
SCW) binned into n=15 bins corresponding to the follow- 
ing energy bands: 3-4, 4-5, 5-6, 6-7, 7-9, 9-11, 11-13, 13-15, 
15-20, 20-30, 30-40, 40-50, 50-70, 70-100, 100-200 keV. 

This PCA showed that only two independent compo- 
nents contribute significantly to the variability of GX 339—4 
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during the state transition. These two components are dis- 
played in Fig.[6l The first component (PCI), responsible for 
75% of the sample variance, consists in a pivoting of the 
spectrum around 6 keV. The second component (PC2) con- 
sists in intensity variation of the hard power-law component 
(almost constant slope) on top of constant soft component. 
It contributes to 21% of the sample variance. The timescale 
of the two variability modes is not firmly determined because 
our observations are not continuous. We can only state that, 
it is comprised between 2 ks, i. e. the minimum duration of 
a pointing, and one month, i. e. duration of the whole ob- 
servation. 

These results are remarkably similar to those obtained 
by Malzac et al. (2006) during an intermediate state of 
Cyg X— 1. These authors found that the spectral variabil- 
ity occurred through two independent modes: the first con- 
sisted in changes in the overall luminosity on time scale of 
hours with almost constant spectrum (fiaring mode); the 
second mode consisted in a pivoting of the spectrum around 
10 keV. The flaring mode was interpreted as fluctuations of 
the heating rate in the corona possibly associated to some 
erratic magnetic activity. The pivoting was associated to a 
mini-state transition. It was interpreted as due to changes 
in the soft cooling photon flux entering the corona as a con- 
sequence of variations of the temperature and luminosity of 
the optically thick disc. 

Unlike Cyg X— 1, in GX 339—4 the pivoting mode is 
more important than the flaring mode, probably because 
Cyg X— 1 was observed during an incomplete spectral tran- 
sition hard-to-intermediate, while GX 339—4 achieved fi- 
nally a softer state. The presence of a pivot at 6 keV dom- 
inating the variability is apparent in the averaged spectra 
shown in Fig. [3] In the framework of the eqpair model, this 
pivoting is clearly associated to the disc getting hotter, with 
the thermal component becoming prominent while the ratio 
Ih/ls was decreasing, leading to a softer and weaker hard 
X-ray spectrum. The interpretation of the pivoting mode in 
both Cyg X— 1 and GX 339—4 as due to variations in the 
soft cooling photon fiux entering the corona is very plausible. 



these results are consistent with the so-called truncated disc 
model (Done et al. 2007). 

In all spectra, including those in of HIMS, we found 
a significant contribution from a non-thermal component. 
This component appears as a high energy excess above the 
pure thermal Comptonisation spectrum. We associated this 
component with the presence of a non-thermal tail in the 
distribution of the Comptonising electron. 

Finally we analysed the variability of the source dur- 
ing the state transition with a PGA. The PCA shows that 
the variability of GX 339—4 can be described as produced 
through two independent modes: 

• spectral evolution with the spectrum pivoting around 6 
keV (PCI) contributing for 75% of the observed variance. 
The pivoting mode may be interpreted as caused by changes 
in the soft cooling photons flux in the hot Comptonising 
plasma associated with an increase of temperature of the 
accretion disc. These changes are well in agreement with 
the variation in the flux of the soft cooling photons found 
by fitting averaged spectra with eqpair. 

• intensity variation of the hard power-law component 
(almost constant slope) on top of constant soft component 
(PC2) contributing at 21%, also consistent with the varia- 
tions of either the compactness ratio Ih/h or the hard X-ray 
luminosity, both found by the spectral fitting. This variabil- 
ity mode could be due to magnetic flares or local instabilities 
in the corona. 
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7 CONCLUSIONS 

We have presented INTEGRAL observations of GX 339-4 
which covered part of the 2004 outburst. We followed the 
source spectral evolution through the Hard-Intermediate 
and Soft-Intermediate states, until the High-Soft state. 

During the HIMS-to-SIMS transition we found that the 
high energy cut-off disappears, in agreement with Belloni et 
al. (2006). Furthermore, the hard X-ray spectrum appears 
softer during the transition states than in the HSS. 

Our detailed spectral analysis of time averaged spectra 
at different stages of the transition confirms that the spec- 
tral transition is driven by changes in the soft cooling photon 
flux in the corona associated with an increase of disc temper- 
atures (leading to dramatic increase of disc luminosity) . The 
measured disc temperature versus luminosity relation sug- 
gests that the internal disc radius decreases. In contrasts, 
the heating rate of the electrons in the corona appears to 
remain nearly constant. Although other models such as dy- 
namic accretion disc corona models cannot be ruled out. 
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Table 1. Log of the INTEGRAL observations. A part the data-set 6 (because of the low exposure time), all periods from 1 to 7 have 
been used for the averaged spectra analysis. The spectral states classification is based on the timing analysis of RXTE data (Belloni et 
al. 2006; Belloni priv. comm.). 



Period 


Rev. 


Start (MJD) 


End (MJD) 


Exp. 


, IBIS 


(ks) 


Exp. JEM-X 


sew 


State 


1 


222 


53226.14 


53227.12 




24.8 




24.3 


6 


HIM 


2 


223 


53228.36 


53229.78 




32.1 




31.8 


8 


HIM 


3 


224.1 


53231.93 


53232.65 




48.5 




48.6 


22 


HIM 


4 


224.2 


53232.68 


53233.13 




30.4 




29.0 


14 


SIM 


5 


226 


53236.61 


53237.77 




33.3 




32.5 


10 


HS 


6 


227-229 


53241.67 


53245.90 




9.07 




6.05 


3 


HS 


7 


232-233 


53255.12 


53259.51 




25.2 




25.4 


13 


HS 



Table 2. Spectral fitting parameters of joint IBIS, JEM-X and SPI (only data-set 1, 2 and 3 for the latter) by using a simple model, 
namely a multicolor-disc black-body plus a power-law with a high energy cut-off (when requested). In periods 1 and 2 the disc temperature 
has been frozen at 0.5 kcV. 



Period 


Rev. 


kTn (kcV) 


r 


Ec (keV) 


Xlidoi) 


1 


222 


(0.5) 


1.4 ±0.1 


71+« 


1.0(243) 


2 


223 


(0.5) 


1.8±0.1 




1.1 (214) 


3 


224.1 


0.70 ± 0.03 


2.1 ±0.1 


iist^J 


1.2(216) 


4 


224.2 


0.69 ±0.03 


2.4 ±0.1 


- 


1.1(173) 


5 


226 


0.64 ±0.02 


1.8 ±0.2 


- 


1.1(121) 


7 


232-233 


0.64 ±0.02 


2.0 ±0.1 


- 


1.0(98) 



Table 3. Best-fit parameters of the joint IBIS/ISGRI, JEM-X and SPI spectra (only rev 222, 223 and 224.1 for the spectrometer). Fits have been performed simultaneously with EQPAIR 
combined with DISKLINE. See text for the parameters description. 



Period 


Rev 


Ih/U 


Inth/lh 


Tes 


kTbfc 


n/2TT 


Ginj 


Ttot 


kTe 


X?(dof) 




Flux xlO 


-9 












[eV] 








[kcV] 




[erg cm 2 c 


''-'] 
























Bol 


bb 


Compt 


1 


222 


4 •^fi+l-02 


75+0-25 
U- '0-0.14 


2 co+0.24 


(300) 


23+0-2'3 
^■^•^-0.13 


„ +0.46 
^•'0-0.64 


2.72 


23.2 


0.99(243) 


11.3 


1.0 


9.4 


2 


223 


9 in+o*" 


10+° 


1 96+°-^^ 


(300) 


0.6 ±0.3 


2 62+0-27 
^•0^-0.37 


2.05 


19.2 


1.17(219) 


12.4 


2.1 


8.9 


3 


224.1 




71+0-18 


n 1S+0-12 

^■l°-0.03 


380+|:^ 


1 n+o-io 

^■U-0.35 


2.9llSi^ 


0.20 


45.0 


1.14(215) 


21.0 


15.4 


5.6 


4 


224.2 


Q ^7+0.03 
'^■-^'-0.02 


'^•o°-o.io 


36+°-^'* 

^■•^"-0.12 


3881^^ 


(1) 


2.861°:?^ 


0.36 


29.4 


1.06(173) 


24.3 


17.8 


6.5 


5 


226 


n 1^+013 


94+0-06 
'J-y^-o.io 


<0.1 


495t^^ 


(1) 


1.9 ±0.3 


0.1 


21.7 


1.15(121) 


13.8 


10.9 


2.9 


7 


232-233 


07+°-''"* 

'^■'J' -0.01 


10+° 

-'-■'^-0.2 


<0.7 


478+66 


(1) 


2.8 ±+0.5 


0.7 


5.5 


1.0(98) 


18.4 


13.8 


4.2 



b 

Co 

O 
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Table 4. Parameters of the SPI and IBIS spectra obtained by fitting data (in 20-500 keV energy range) with thermal-Comptonisation 
models and thermal-Comptonisation plus power-law. COMPPS and COMPTT have been used as therinal-Coinptonisation models. Comp- 
tonisation temperature (fcT) and Thomson depth (t) are free parameters; F is the power-law photon index. In the COMPPS model the 
black body temperature of the soft seed photons was fixed at 0.1 keV in all fits. Reflection component (parameter Q/2n) is only present 
in COMPPS. 



Rev 



Model kT (kcV) 



tau-y 



n/2TT 



xl(dof) F-test p 



222 



COMPTT 


27 ±2 


1.7±0.1 


- 


- 


COMPTT+PO 


21 ±3 


9 1 +0.6 


- 


- 


COMPPS 


42tt 


- 


3.1 ±0.2 


0.6 ±0.3 


COMPPS+PO 


27tl!> 


_ 


tr 2+4.8 


n 4+0-6 



1.8! 



2.0Z 



1.2(111) 

1.1(109) 9.2 X 10" 

1.08(104) 

1.05(102) 0.1 



223 



COMPTT 
COMPTT+PO 

COMPPS 
COMPPS+PO 



561 
17l 



20 



hio 

-6 

+31 



0.6 + 0.3 



2.6 



+6.1 



1.1 + 0.1 



0.7: 

0.2 



+2.4 



2.2 



2.08^ 



+0.8 



.1(104) 


- 


.0(102) 


3.4 X 10"^ 


.2(103) 


- 


.1(101) 


6.6 X 10-2 



